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ABSTRACT 

To understand the influence of additional wide stellar companions on planet forma- 
tion, it is necessary to determine the fraction of multiple stellar systems amongst the 
known extrasolar planet population. 

We target recently discovered radial velocity exoplanetary systems observable from 
the northern hemisphere and with sufficiently high proper motion to detect stellar 
companions via direct imaging. We utilize the Calar Alto 2.2 m telescope in combina- 
tion with its lucky imaging camera AstraLux. 

71 planet host stars have been observed so far, yielding one new low-mass 
(0.239 ± 0.022 M0) stellar companion, 4.5arcsec (227 AU of projected separation) 
northeast of the planet host star HD 185269, detected via astrometry with AstraLux. 
We also present follow-up astrometry on three previously discovered stellar compan- 
ions, showing for the first time common proper motion of the 0.5 arcsec companion to 
HD 126614. Additionally, we determined the achieved detection limits for all targets, 
which allows us to characterize the detection space of possible further companions of 
these stars. 
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1 INTRODUCTION 

Recent radial velocity (RV) and transiting exoplanet sur- 
veys have resulted in the discovery of a large number of 
new exoplanetary systems. Most stars are born in binary o r 
multiple systems l|Mathieu et al.|[2000l . lDuchene et aLll2007l ). 
and indeed we observe that a reasonable fraction of the 
stars in the Galaxy are multiple l|Ladall2006t ). The effects of 
additional stellar companions on the planet formation pro- 
cess are therefore of high interest to constrain and calibrate 
planet formation theories. 

A number of studies have been co nducted on this s ubject 
in the past, such as de scribed in iMugrauer et all (l2007l) 



Daemgen et al.l (|2009|) and m ost recently in Roll et all 
l|201ll ) and lChauvin et all (|201ll '). As a result of these stud- 



ies, 43 multiple stellar systems hosting exoplanets are known 
to date, suggesting that about 17 % of all kno wn exoplanets 
reside in such systems (|Mugrauer et al.ll2009l ). 
In this paper we present the results of our ongoing multi- 
plicity study of the most recently (since 2008) discovered RV 
exoplanet host stars. We utilized the Cal ar Alto 2.2 m tele - 
scope and the AstraLux instrument ( Hormuth et al.ll2008l ) 
to achieve imaging data with higher Strehl ratio compared 
to simple seeing-limited imaging. 

In the following section we characterize our sample. In sec- 
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tion [3] we give a brief introduction to the observation tech- 
nique as well as the instrument used and describe the re- 
duction and astrometric calibration of our data. In section 
[4] we show all confirmed or rejected companion candidates 
with the associated proper motion analysis. We also present 
detection limits for all studied systems. Finally in section [S] 
we summarize our findings. 



2 TARGET SAMPLE 

Our sample consists of stars with RV planet candidates dis- 
covered between 2008 and 2011. They are all observable 
from the northern hemisphere with declinations down to 
-22° and a relatively even distribution in right ascension. 
Since we want to test companionship of detected companion 
candidates via astrometry, we chose targets which are close 
enough to show a sufficient proper motion (e.g. > 47mas/yr 
or approximately one Astralux pixel per year) to do a com- 
mon proper motion analysis approximately one year after 
the first epoch images were taken. The distance distribution 
of our targets is shown in Fig. \T\ On average our targets 
are 56.6 pc away and show a proper motion of 180.4 mas/yr 
equivalent to 3.8 pix/yr on the Astralux instrument. 
The RV planet search technique favors main-sequence stars, 
therefore the average age of our target stars is 4.6 Gyr (data 
as listed in exoplanet. eu bv lSchneider et al]l201ll '). Spectral 
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Figure 1. Histogra m of distances of all ou r target stars (as listed 
in cxoplanet.eu by ISchneider et al.ll2011^ . The majority being 
within 100 pc. 



Figure 3. Minimum mass versus separation for all planet can- 
didates orbi ting hosts included in our sample (as listed in cxo- 
planet.eu bv lSchneider et al.ll20lil ). For comparison we also plot 
minimum mass versus separation of all RV planet candidates dis- 
covered so far. 
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Figure 2. Metallicity distribution for all targets in our sample. 
For comparison we plot the metallicity distribution of all RV 
planet candi date host stars discov ered so far (as listed in cxo- 
planet.eu bv lSchneider et al.|[201l[) . 



Figure 4. Eccentricity distribution of planet candidates orbiting 
hosts included in our sample. For comparison wc also plot the 
eccentricity distribution of all RV planet candidates di scovered 
so far (as listed in cxoplanet.eu bv lSchneider et al.ll20lil ). 



types range from late F to early M, with the majority being 
G and K type. In Fig.[2]we show the metallicity distribution 
of our sample compared to the whole RV planet candidate 
host population. The average metallicity of our sample is 
slightly higher than solar metallicity at 0.09. The metallic- 
ity distribution of our sample does not differ significantly 
from the whole population. 

In Fig. [3] we present the minimum masses of detected RV 
planet candidates versus their semi-major axis for our sam- 
ple and the whole population. The majority of the detected 
planet candidates are located within 2AU of their hosts. 
The average minimum mass of our sample is 2.8Mjup, as 
compared to 2.4Mjup for the whole RV planet candidate 
population. Given the standard deviation of 3.3 M,j„p and 
3.8Mjup respectively this is not significant. 
The eccentricity distribution of RV planets around our sam- 
ple stars can be seen in Fig. (4] We compared the eccentricity 
distribution of our sample to that of all RV planet candidates 
and found no significant differences. 

We conclude that the general properties of our sample re- 
semble those of the whole population. 



3 OBSERVATIONS OUTLINE AND DATA 
CALIBRATION 

3.1 Lucky imaging technique 

To achieve high spatial resolution as well as high sensitivity, 
one must consider the effects of turbulent atmosphere. Given 
the Rayleigh criterion for the minimum resolvable angle of 
— 1.22- A-D~^, the Calar Alto 2.2m telescope has a theo- 
retical resolution of 89 mas at a wavelength of 776 nm, which 
is the central wavelength of the SDSS i' filter which we are 
using. The median seeing at the Calar Alto site, however, 
is 0.9 arcsec l|Sanchez et al.ll2007l ). which is about 10 times 
larger. 

The atmospheric conditions in the optical regime are sub- 
ject to rapid variations with a timescale that can be approx- 
imated by the speck le coherence time Te ~ 0.36 ■ ro ■ Ai/~^ 
(|Roddier et al.lll982l ). where ro is the Fried-parameter and 
/S.V is the wind speed dispersion in the atmosphere. For a 
typical wind speed dispersion of A.v ~ 10m-s~^ at the Calar 
Alto site a nd a V-band seeing of 0.7 arcsec, Te is in the order 
of 100 ms (iHormuth. Diploma Th"esi3 120071 '). Within the 



© 2011 RAS, MNRAS 000, [THM] 



A lucky imaging multiplicity study of exoplanet host stars 3 



resulting speckle pattern will remain fixed, whereas longer 
integration times would lead to an averaged and therefore 
"smeared out" speckle pattern. 

The lucky imaging approach consists of taking several thou- 
sand short images with integration times shorter than Te , to 
sample the speckle variations during the observation win- 
dow. We then only choose the so called "lucky shots" with a 
very high Strehl ratio in one of the speckles, to shift and add, 
resulting in a final image with the highest possible strehl ra- 
tio and therefore highest possible angular resolution. 
For an in- depth intro duction to the lucky imaging technique 
please see iLaw et al.l ((2OOQ). 



3.2 Observations and data reduction 

All observations were carried out with the Calar Alto 2.2 m 
telescope in combination with the AstraLux instrument. 
This instrument consists of a back-illuminated, electron- 
multiplying, frame transfer CCD, which is well suited for 
the lucky imaging observation t echnique. For a detailed de- 
scription of the instrument, see iHormuth et all (|2008l ). 
Since we wanted to detect low-mass stellar or even brown 
dwarf companio ns of our target s, we ch ose the SDSS i' filter 
as described in iFukugita et al.l (|l996l l for all observations. 
This is the best choice taking into account the sensitivity 
of the detector and the brightness of the detectable com- 
panions, as well as the variability of the atmosphere. For 
calibration we took dome- and sky-flats at the beginning 
and/or end of each observation night. 

For all science targets we chose 29.54 ms of exposure time 
per frame, which is well below the typical speckle coherence 
time at the Calar Alto site. Shorter exposure times would 
have led to signiflcantly growing overheads since we would 
have needed to switch off the frame transfer mode of the 
instrument. Longer exposure times would have resulted in a 
less localized and more "smeared out" speckle pattern. 
The electron-multiplying gain was individually adjusted for 
each target to obtain the maximum signal while still oper- 
ating in the linear regime of the detector. We then took 100 
dark frames with closed shutter and with the respective elec- 
tron multiplying settings before the start of each acquisition 
sequence. 

Depending on weather conditions and time constraints we 
took between 10000 and 80000 frames per science target, 
which corresponds to a total integration time of 0.25 min 
and 1.97 min respectively for a typical frame selection rate 
of 5%. The obtained AstraLux images are all reduced and 
processed with our own ESO-MIDAS (Munich Image Data 
Analysis System) software for the reduction of imaging data 
taken with the lucky imaging technique. The individual As- 
traLux images are dark-subtracted with the median (for op- 
timal suppression of cosmic ray artifacts) of several short in- 
tegrated dark-frames and flat-fielded. For the lucky imaging 
data processing, the reached Strehl-Ratio in all individual 
images is determined at first using the bright planet host 
star as Strehl-probe. The frames are then ranked accord- 
ing to their Strehl-Ratio and only the best 10, 5, and 1 % 
of all frames, i.e. only the frames with the highest Strehl- 
Ratio, are then selected and combined using the shift -l-add 
technique. Thereby, for image registration, the position of 
the brightest pixel in the speckle pattern of the planet host 



Table 1. Hipparcos astrometry of epoch 1991.25 of our calibra- 
tion binaries 



Binary 


Separation [arcsec] 


Position Angle [°] 


HIP 59585 


18.677 ± 0.032 


190.83 ±0.10 


HIP 65205 


15.041 ± 0.015 


219.44 ± 0.06 


HIP 67099 


18.751 ± 0.028 


343.33 ± 0.09 


HIP 72508 


15.248 ± 0.014 


91.36 ±0.05 


HIP 80953 


16.356 ± 0.021 


195.69 ±0.07 



Table 2. Astrometric calibration of all observation epochs. We 
list the pixel-scale (PS) and the position angle (PA) of the y-axis 
for all observation epochs. 



Epoch 


PS [mas 


/pix] 


PA of y-axis [°] 


23/04/2008 


47.278 


± 


0.101 


0.17 ±0.23 


11/07/2008 


47.180 


± 


0.078 


0.44 ±0.15 


16/01/2009 


47.237 


± 


0.114 


359.39 ± 0.08 


07/09/2009 


47.220 


± 


0.100 


0.43 ±0.13 


23/02/2010 




_i 




_i 


14/07/2010 


47.243 


± 


0.049 


358.49 ± 0.22 


14/01/2011 


47.365 


± 


0.135 


358.07 ±0.12 


27/07/2011 


47.160 


± 


0.066 


358.37 ±0.23 


^ due to bad weather conditions no astrometric cal- 



ibrators could be observed 



star is determined in all frames, which are then shifted and 
averaged. 

3.3 Astrometric calibration 

Since we want to precisely measure separation and position 
angle of our companion candidates with respect to the pri- 
mary stars, and the AstraLux instrument is not permanently 
mounted to the Calar Alto 2.2m telescope, a careful astro- 
metric calibration of the detector in each observation epoch 
is necessary. For this purpose we selected a sample of five 
wide binaries, for which precise Hipparcos measurements as 
well as several additional observation epochs are available. 
We show our sample along with the Hipparcos astrometry 
in Table [T] Furthermore, we observed the center of the glob- 
ular cluster M15 whenever possible, for which precise HST 
astrometry is available. 

For each epoch we have a minimum of three calibrators 
which were observed in the same night as the science tar- 
gets. For all calibrators we used the Hipparcos measurement 
of epoch 1991.25 as reference point and then linearly fitted 
the slow orbital motion by using all data points available in 
the Washington Double Star (WDS) catalog. 
The pixel scales and position angles of the calibration images 
for each individual calibrator were calculated using ESO- 
MIDAS for the position measurements of the binary com- 
ponents, and GAIA (Graphical Astronomy and Image Anal- 
ysis Tool) for the respective measurements of the M15 clus- 
ter compo nents (using SExtractor: S oftware for source ex- 
traction bv lBertin and Arnouts|[l996l ). The average of these 
calculations were used as final calibration for the respective 
epoch to cancel out systematic errors due to residual orbital 
motion of our binary stars. The final astrometric calibrations 
for all observation epochs are listed in Table (2] 



© 2011 RAS, MNRAS 000,[lHT4] 



4 C. Ginski et al. 



3.4 PSF subtraction and astrometric 
measurements 

Since we are searciiing for faint low-mass stellar compan- 
ions to bright stars, we subtract the primary stars' PSF in 
order to detect the companions with highest possible signal- 
to-noise. Due to the random nature of the lucky imaging 
approach, the PSFs in the final reduced images may vary 
significantly from star to star, for example as a function of 
the average zenith distance of the observed objects. Hence 
it is not feasible to subtract a PSF standard from all the 
images. Furthermore, the PSFs of our target stars are not 
necessarily symmetrical in nature, rendering simple rotation 
subtraction techniques ineffective. 

Given these limitations, we found the best approach to be an 
unsharp mask filter. We fold the images with a 2-dimensional 
Gaussian of size n to "blur" them, and then subtract the 
"blurred" images from the originals, thus eliminating all low 
spatial frequencies from the images. We tested sizes n be- 
tween 5 and 20 pixels for different images, with n = 10 
generally yielding the best results in terms of improving the 
signal-to-noise of the faint companion candidates. 
We again used ESO-MIDAS for all astrometric measure- 
ments, first measuring the primary star's position in the 
original images, then applying the PSF subtraction as de- 
scribed, before measuring the companion candidate's posi- 
tion. We always measured both positions multiple times, 
averaging the results to avoid statistical errors. 



4 RESULTS 

To date we observed 71 planet host-stars in our project. 
In seven cases we detected faint companion candidates, of 
which three were already known. In the following para- 
graphs, we show the astrometric measurements for all the 
companion candidates, as well as the proper motion analysis, 
to determine which of the objects are physically associated 
to the target stars. 

4.1 Imaging of known companions 

In the course of our study we imaged three already known 
companions to exoplanet hosts, to test the instrument capa- 
bilities and to provide additional astrometric measurements 
for future orbit solutions. The final reduced images can be 
seen in Fig. [5] and our astrometric measurements are listed 
in Table El 

The stellar companion to r Boo resides at an angular sep- 
aration 2.181 ± 0.018 arcsec, which corresp o nds t o a pro- 
jected separation of 34 AU. [Patience et al.l l|2002l ') stated, 
that T Boo A and B form an eccentric binary system with a 
semi-major axis of « 225 AU, and masses of « 1.3 M© and 
« 0.4 respectively. 

HD 176051 A and B form a close binary system of only 
1.139 ± 0.009 arcsec separation and an orbital period of 
61.2 yr . The stellar masses are 1 . 07 Mp, and 0.71 M© respec - 
tively (|Muterspaugh et al.llioToI 'l. lMuterspaugh et al.l (|2010l ') 
found an astrometric companion of ~ 1.5 Mj, to be orbiting 
one of the components. 

The stella r companion to HD 1 26614 was only detected re- 
cently by [Howard et al.l (|2010l ) at the Palomar Observa- 



Table 3. Astrometric measurements of all known eompanions 



Primary 


Epoch 


Separation 


Position Angle 






[arcsec] 


[°] 


T Boo 


23/04/2008 


2.181 ±0.011 


46.26 ±0.42 


HD 176051 


27/07/2011 


1.139 ±0.005 


252.12 ± 0.40 


HD 126614 


14/01/2011 


0.499 ± 0.067 


60.70 ±5.60 



tory, using adaptive optics imaging. It is separated from 
the primary only by 0.490 ± 0.051 arcsec, which corresponds 
to a projected separation of 35.6 AU. The mass of the 
companion was inferred from the JHK photometry, to be 
0.324 ± 0.004 Mq. Its companionship was so far only con- 
firmed by photometry and radial velocity trend. We are not 
aware of any further observations of this companion, es- 
pecially regarding whether or not it is actually co-moving 
with the primary. Hence we present a proper motion analy- 
sis in section 14.31 using the measurements bv lHoward et al.] 
(| 2010! ) and our own AstraLux data point, to confirm that 
HD 126614 A and B form indeed a common proper motion 
pair. 

4.2 New companion candidates 

In the course of our study we found new faint compan- 
ion candidates around the planet host stars HD 185269, 
HD 183263, HD 187123 and HD 13931. In Fig. [6] we show 
the epochs with the highest signal-to-noise. 

In Table |4] we show all the astrometric measurements 
obtained in the different observation epochs for all compan- 
ion candidates. Given the distances of the primary stars, the 
detected objects would have projected separations ranging 
from 179 AU (close candidate of HD 187123) up to 781 AU 
(far candidate of HD 187123). 

4.3 Proper motion analysis 

We performed a proper motion analysis for all companion 
candidates listed in Table |4] to distinguish between true 
physical companions and distant background objects. Hence 
we computed how separation and position angle would 
change from the first to the last observation epoch if the 
detected objects are indeed only in the background, and 
are therefore most likely standing still, also considering 
the effect of the annual parallax due to the earth's motion 
around the sun. We then compared these expected positions 
with the actual positions measured in our images. The 
resulting diagrams for all new companion candidates are 
listed in Fig. [T] Even if the detected objects are physically 
associated to the primaries, they might show a slightly 
different proper motion, which can be explained by orbital 
motion around the primary star. To account for this, we 
also added to the analysis the maximum possible change 
of separation and position angle due to orbital motion at 
the projected distance of our first epoch measurement. We 
assumed an edge-on circular orbit for changes in separation 
and a face-on circular orbit for changes in position angle 
respectively. 

HD 13931 
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(a) T Boo taken in epoch 23/04/2008 (b) HD 176051 taken in cpoeh 27/07/2011 (c) HD 126614 taken in epoeh 14/01/2011 

Figure 5. Follow-up images of already known companions to exoplanet host stars taken in SDSS i' filter during our ongoing campaign. 
All images have been PSF subtracted. North is always up and east to the left. The companions are marked by arrows. 



Table 4. Astrometric measurements of all companion candidates 



Primary Epoch # CC Separation [arcsec] Position Angle [°] 



HD 13931 


14/01/2011 


1 


9.972 ± 0.030 


19.50 ±0.16 




27/07/2011 


1 


10.003 ± 0.015 


19.16 ±0.25 


HD 185269 


07/09/2009 


1 


4.522 ± 0.014 


7.91 ±0.23 




14/07/2010 


1 


4.520 ±0.010 


8.17 ±0.30 




27/07/2011 


1 


4.511 ± 0.008 


8.44 ±0.28 


HD 183263 


23/04/2008 


1 


14.190 ± 0.032 


47.83 ±0.26 






2 


3.295 ± 0.015 


349.55 ±0.40 




11/07/2008 


1 


14.183 ±0.028 


47.92 ±0.19 






2 


3.245 ±0.0171 


349.72 ±0.3li 




07/09/2009 


1 


14.244 ± 0.033 


47.93 ±0.17 






2 


3.315 ± 0.031 


350.75 ±0.56 


HD 187123 


11/07/2008 


1 


2.926 ±0.011 


48.11 ±0.29 






2 


11.560 ± 0.022 


263.12 ±0.19 




07/09/2009 


1 


2.917 ±0.013 


43.92 ± 0.29 






2 


11.696 ± 0.028 


263.92 ±0.17 



^ due to the spurious detection of the companion candidate this data point seems 
dubious 



HD 13931 was observed in January and July 2011. The 



associated proper motion plot can be found in Fig. 7(a 



The data point taken in July 2011 is, within its uncertainty, 
consistent with the background and common proper motion 
hypothesis. The common proper motion hypothesis can be 
rejected only on the 0.83 and 1.09 cr level for separation and 
position angle respectively, hence it is not possible to draw 
a final conclusion from our AstraLux measurements. 
This companion candidate was resolved in 2MASS (Two 
Micron All Sky Survey) J, H and Ks-band images taken on 
October 28th 1998. We show the corresponding Ks-band 
image in Fig. [S] Our measurements for this epoch yield 
a separation of 8.13 ± 0.15 arcsec and a position angle 
of 31.14° ± 0.75°. Using this new data point, we created 
a proper motion plot as shown in Fig. [9] Although the 
2MASS data points are not consistent with the background 
hypothesis within Ict, we can still reject common proper 
motion with 10.38 cr for separation and 11.29 a for position 
angle. The inconsistency with the background hypothesis is 
most likely caused by a slow proper motion of the detected 
object. We conclude that our companion candidate around 




Figure 8. 2MASS Ks-Band Image of HD 13931 in epoch 
28/10/1998. The AstraLux field of view is indicated by the white 
rectangle. 



HD 13931 is only a background object. 
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(a) Image of HD 13931 taken in epoch 14/01/2011 



(b) Imago of IID 185269 taken in epoeli 07/09/2009 




(c) Image of HD 183263 taken in epoch 23/04/2008 



(d) Image of HD 187123 taken in epoch 07/09/2009 



Figure 6. Images of our newly discovered companion candidates taken in SDSS i' filter. Images are 24.2 X 24.2arcsec, north is always 
up and east to the left. In all cases the primary stars PSF was subtracted by an unsharp mask filter. Companion candidates are marked 
with circles and arrows. The numbers assigned to multiple companion candidates match those used in Table |4] 
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(c) Companion candidate 1 to HD 183263 
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(b) Companion candidate to HD 185269 
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(d) Companion candidate 2 to HD 183263 
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(e) Companion candidate 1 to HD 187123 (f ) Companion candidate 2 to HD 187123 

Figure 7. Proper motion diagrams for all newly discovered companion candidates. Upper diagrams are always for the position angle 
and lower diagrams for the separation, both plotted over time in years. The dashed lines mark the area for maximum possible orbital 
motion in case of bound objects and circular orbits. The solid lines mark the area for non-moving background objects. The wobble in 
the background hypothesis is introduced by the annual parallax. 
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Figure 9. Proper motion diagram for the companion candidate to 
HD 13931 including a 2MASS measurement of epoch 28/10/1998. 
The areas for maximum orbital motion and background hypoth- 
esis are plotted backwards from our latest measurement in epoch 



27/07/2011. 



HD 185269 



In Fig. 7(b) we show the proper motion plot for the 
companion candidate next to HD 185269. Due to the 
overlap of background and common proper motion area 
in the position angle plot, our angular measurements 
would be consistent with both hypothesis. However, our 
measurements of the separation are only consistent with an 
object of common proper motion with the primary star. We 
can reject the background hypothesis with 9.9 cr, and are 
therefore concluding that our object is indeed a previously 
unknown stellar companion to HD 185269. HD 185269 B 
has a projected separation of 227 AU . Further follow-up 
observations will be executed to determine the orbital 
motion of the B component. 

Using IRAF (Image Reduction and Analysis Facility) 
standard aperture photometry, we calculate an average 
magnitude difference of A/ = 7.2 ± 0.2 mag throughout 
all our obser vation epochs. Given the age of the primary 
of 4.2 Gyr Jjohnson et all 120061') and the parallax of 



21.11 ± 0.74 mas ( van LeeuwenI 20071 ), we find a mass o f 



0.239 ± 0.022 Mq, using the models bv lBaraffe et all l|l998f) . 



HD 183263 



Fig. 7(c) and 7(d) show the proper motion plots for the two 
companion candidates discovered around HD 183263. 
Due to the comparatively small proper motion of only 
38.09 mas-yr~^ (~ O.SAstraLux pix), a final conclusion 
for the first companion candidate is difficult with our 
epoch difference of 1.25 yr. The position angle gives no 
useful information, since the background and common 
proper motion area are completely overlapping, and we 
can only reject the common proper motion hypothesis with 
1.06 (T, in the case of the separation. We could, however, 
identify our first companion candidate once more in the 
2MASS catalogue. We show the corresponding image in 
Ks-band in Fig. 1101 The image was taken on August 5th 
1999, giving us 10 yr of epoch difference. The 2MASS-PSC 



Figure 10. 2MASS Ks-Band Image of HD 183263 in epoch 
05/08/1999. The AstraLux field of view is indicated by the white 
rectangle. 




1998.0 2000.0 2002.0 2004.0 2006.0 2008.0 2010.0 
time [year] 



Figure 11. Proper motion diagram for the first (far) compan- 
ion candidate to HD 183263 including a 2MASS measurement of 
epoch 08/1999. The areas for maximum orbital motion and back- 
ground hypothesis are plotted backwards from our latest mea- 
surement in epoch 07/09/2009. 



(2MASS - Point Source Catalogue, ICutri et al.1 |2003| ) 
measurements yield a separation of 13.8 ± O.larcsec and 
a position angle of 48.2° ± 0.4°. Using this information, 
we created a new proper motion plot incorporating the 
2MASS epoch which is shown in Fig. 1111 Even though 
the error bars are significantly larger than the ones of our 
AstraLux measurements, the new proper motion analysis 
in separation clearly indicates that the first companion 
candidate around HD 183263 is a background object, with 
a significance of 3.3 cr. 

The second (closer) companion candidate was also detected 
in all three observation epochs. In July 2008 it was only 
barely visible, making these data points somewhat dubious. 
While the separation measurements are in general still 
consistent with a co-moving object as well as a background 
object, the development of the position angle indicates 
(with 1.4 a corresponding to 83.8% probability) that the 
second companion candidate is a background object as well. 
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2009.5 2010.0 2010.5 

time [year] 



Figure 12. Proper motion diagrams for tiic stellar companion to 
HD 126614 



HD 187123 



The proper motion plots for the two companion can- 



didates around HD 187123 are shown in Fig. 7(e) and 7(f) 



Both companion candidates show as fully consistent with 
the background hypothesis with an average significance 
level of 3.3 a. We conclude that both objects are indeed 
distant, non-moving background objects. HD 187123 is a 
single star within the given detection limits described in 
section 14.41 

HD 126614 

We also performed a proper motion analysis for the 
companion to HD 126614, using the measurements by 
[Howard et al.l (|2010l ) as first epoch, and our own AstraLux 
measurement as second epoch. The resulting plot is shown 
in Fig. [H 

Our measured separation of 499 ± 67 mas and position 
angle of 60.7° ± 5.6° are consistent with the the astrometry 
of the companion in the dis covery epoch (489.0 ± 1.9 mas 
and 56.1° ±0.3°, as given bv [Howard eralll2010l ) . However, 
due to the very low declination of —05° 10' 4", we had to ob- 
serve HD 126614 at high air mass, which led to an extended 
PSF in north-south direction. This made it very difficult to 
determine the position angle precisely, which is also evident 
in the large error bars of the measurement. Still the posi- 
tion angle is, within its error bars, consistent with common 
proper motion. 

Despite the uncertainties in the position angle, we can still 
reject the background hypothesis with 1.7 cr and 5.3 ct for 
position angle and separation respectively. We can therefore 
conclude that HD 126614 B is indeed physically associated 
with A. 

Due to the small separation between primary and com- 
panion, aperture photometry with IRAF is not viable. 
We could, however, infer the I-band magnitude differ- 
ence of companion and primary from the ratio of the 
peak counts. We c alculate AJ = 5.59 ± 0.15 mag. Given 
the log(age) = 9.1 {Howard ct al. 2010*) and parallax of 
14.63 ± 1.15 mass (|van Leeuwen .2007 ) , we find a mass o f 
0.307 ± 0.033 M0, using the models bv lBaraffe eTsI\ l|l998l ). 



This is consis t ent w ith the mass of 0.324 ± 0.004 found by 
[Howard et al] l|2010l) . 



4.4 Detection limits and non-detections 

In Table [6] and Table [7] we list all the systems that were 
observed in our survey to date, with and without additional 
stellar components detected, respectively. We also list the 
epochs in which these systems have been observed, since 
our detection limits vary with the weather and seeing con- 
ditions. Additionally we present limits for the lowest mass 
objects detectable at separations of 0.5, 1, 2 and 5arcsec, at 
a signal-to-noise ratio of 5, in the PSF subtracted images. 
We calculated the uncertainties of these mass limits consid- 
ering the uncertainties in the measured parallaxes for each 
system, as well as a maximum photometric error of 0.2 mag 
in our Astralux images, along with an average photometric 
error of 0.05 mag in the primary stars' I-band magnitudes. 
In Fig. [13] we show the average dynamic range plots for each 
of our observation epochs. We determined the average dy- 
namic range by measuring the noise levels in all reduced 
images of a given epoch with a frame selection rate of 5 %, 
assuming a signal-to-noise ratio of 5. The primary stars' PSF 
was always subtracted to get the best possible results. 
In all but two epochs we reach a magnitude difference of 
A7 = 10 mag outside of 2arcsec from the primary star. 
Given the average absolute magnitudes and ages of the tar- 
get stars in each epoch as listed in Table [S] (calculated from 
individual values as listed in Simbad and exoplane t . eu re - 
spectively), and using the models bv [Baraffe et al.[ l|l998l ). 
we would have easily been able to detect all wide stellar 
companions. Furthermore, we still probe the low mass stel- 
lar regime (M < 0.15Mq) down to separations of larcsec 
from the primary star. 

The observation epochs in January 2009 and February 2010 
suffered from bad weather conditions, and while we could 
complete most of our observation program in January 2009, 
we could only observe for approximately one hour through 
the cloud cover in February 2010. This is reflected in the 



dynamic range graphs 13(c) and 13(e) which show a no- 



ticeable decline in sensitivity in comparison with the other 
observation epochs. On average we are still able to detect 
stellar companions down to 0.15 A^q outside of 1.2arcsec 
and 3.1arcsec respectively. 

If we average all minimum masses for the respective separa- 
tions as derived in Table [T] we conservatively get detection 
limits of 0.34 Mq at 0.5 arcsec as well as 0.22 Mq at 1 arcsec 
and 0.15 Mq outside of 2 arcsec. Given our field of view, we 
can detect such low-mass stellar companions up to separa- 
tions of 12 arcsec. 



5 CONCLUSIONS 

In our ongoing study we observed 71 planet host stars to 
date. Of these 71 systems, 3 were already known to be mul- 
tiple, for which we present follow up astrometry. Thereby, 
we show for the first time conclusively, that the compan- 
ion to HD 126614 is indeed physically associated with the 
primary. As the planet host star HD 126614 also exhibits 
a further companion at a wider separation (HD 216614 C, 
sep= 41. 914±0. 110 arcsec (3043 AU),PA= 299.36° ±0.14° 
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(a) average dynamic range in epoch 23/04/2008 (b) average dynamic range in epoch 11/07/2008 





(c) average dynamic range in epoch 16/01/2009 (d) average dynamic range in epoch 07/09/2009 





(e) average dynamic range in epoch 23/02/2010 (f) average dynamic range in epoch 14/07/2010 






(g) average dynamic range in epoch 14/01/2011 (h) average dynamic range in epoch 27/07/2011 



Figure 13. Average dynamic range plots for a signal-to-noise of 5 for all our observation ep ochs. The das h ed lin es represent the detectable 
minimum mass objects at a given separation from the primary star, using the model by iBaraffe et all l ll99St) . The primary stars' PSF 
was always subtracted by unsharp mask filtering before determining the dynamic range. 
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Table 6. List of all stars with companions and companion candidates detected. We give the numbers of frames observed, each frame 
with an e xposure time of 29.5 4 ms. We also give an upper limit for the minimum mass detectable at 0.5, 1, 2 and 5arcsec using the 
models bv lBaraffe et al] lll998l ) and assuming an average age for each epoch as listed in Table [5] 



Star 


Epoch 


# of frames 


M [Mq] at 0.5" 


M [Mq] at 1" 


M [Mq] at 2" 


M [Mq] at 5" 


HD13931 


14/01/2011 


50000 


0.173 ± 0.014 


0.097 ± 0.003 


0.086 ± 0.001 


0.085 ± 0.001 




27/07/2011 


50000 


0.365 ± 0.036 


0.142 ± 0.010 


0.105 it 0.004 


0.100 ± 0.004 


HD183263 


23/04/2008 


50000 


0.194 ± 0.020 


0.099 ± 0.003 


0.086 ± 0.001 


0.085 ± 0.001 




11/07/2008 


80000 


0.180 ± 0.016 


0.144 ± 0.011 


0.099 ± 0.003 


0.088 ± 0.002 




07/09/2009 


80000 


0.178 ± 0.015 


0.096 ± 0.003 


0.087 ± 0.001 


0.086 ± 0.002 


HD187123 


11/07/2008 


77800 


0.149 ± 0.010 


0.103 ± 0.004 


0.086 ± 0.001 


0.083 ± 0.001 




07/09/2009 


80000 


0.160 ± 0.012 


0.091 ± 0.002 


0.084 ± 0.001 


0.084 ± 0.001 


HD185269 


07/09/2009 


50000 


0.307 ± 0.032 


0.166 ± 0.013 


0.107 ± 0.005 


0.096 ± 0.003 




14/07/2010 


60000 


0.286 ± 0.027 


0.125 ± 0.007 


0.093 ± 0.003 


0.091 ± 0.002 




27/07/2011 


50000 


0.239 ± 0.021 


0.113 ± 0.006 


0.091 ± 0.002 


0.089 ± 0.002 


T Boo 


23/04/2008 


50000 


0.208 ± 0.019 


0.104 ± 0.004 


0.088 ± 0.001 


0.088 ± 0.001 


HD176051 


27/07/2011 


50000 


0.149 ± 0.010 


0.095 ± 0.003 


0.084 ± 0.001 


0.084 ± 0.001 


HD126614 


14/01/2011 


50000 


0.179 ± 0.018 


0.103 ± 0.005 


0.090 ± 0.002 


0.088 ± 0.002 



Table 5. Average absolute magnitudes and ages of all target stars 
in an observation epoch 



Epoch MI [mag] Age [Gyr] 



23/04/2008 


4.32 


4.9 


11/07/2008 


3.43 


6.0 


16/01/2009 


2.93 


4.6 


07/09/2009 


2.64 


6.2 


23/02/2010 


3.72 


2.2 


14/07/2010 


4.14 


7.4 


14/01/2011 


3.67 


5.4 


27/07/2011 


3.20 


6.5 



at 2MASS Epoch May 3th 2000), the HD 126614 system is 
actuaUy a hierarchical triple system, the only one presently 
known in which the planet host star exhibits a close stellar 
companion. All other known planet host triples are com- 
posed of the planet host star and a binary-companion at 
wider separation. 

We also discovered one new low-mass (0.239 ± 0.022 M©) 
stellar companion to the star HD 185269, with a separation 
of 4.511 ± 0.013 arcsec at a position angle of 8.44° ± 0.30°. 
This corresponds to a projected distance of 227 AU. 
HD185269A harbors a "Hot Jupiter" {M sini = 0.94 Mj) 
with an orbital peri od of 6.8 d and a se mi-major axis of 
0.077 AU, detected bv ljohnson et all (j2006l ). They state that 
its orbital eccentricity of 0.3 is large in comparison with 
other planets found within 1 AU of their host stars. This 
could poss ibly be explained by the Kozai mechanism as de- 
scribed bv lTakeda and Rasiol l|2005l ). We used the formula 

P - 2nJ . . fH^^) (^) ' (1 - elf' 

V G{mo + mi) V m2 J \aiJ ^ ' 

to ca lculate the period of the Kozai oscillations (iFord et al.l 
I2OOOI ). where the indices 0, 1 and 2 represent the host star, 
planetary and stellar companions respectively. Assuming an 
eccentricity of 0.5 for the stellar compan ion (as suggested by 
statistical analysis in ISoderhielnj|l999l ). and a semi-major 
axis equal to the projected separation of 227 AU, we get a 
period of 1.7 Gyr. Given the age of HD 185269 A of 4.2 Gyr, 
this period is short enough so that the Kozai effect might 



have altered the planetary companion's eccentricity and in- 
clination. 

For the remaining 66 target stars in our sample we can 
on average exclude all low-mass stellar companions of 
M > 0.15 Mq down to 2 arcsec around the primary. For 
many stars in our sample we get significantly deeper, en- 
abling us to exclude all low-mass stellar companions outside 
of 2 arcsec and up to 12 arcsec. To further summarize our 
detection limits, in Fig. [14] we show the average dynamic 
range plot of all observations, with all detected stellar com- 
panions marked. 

68 stars of our sample had an unknown multiplicity status 
before our observations. Of those, only one proved to be mul- 
tiple within the detection limits of our survey. This yields 
a multiplicity rate of only 1.5 % for our sample. Given that 
the multiplicity rate of the known exoplanet host population 
is about 17%, stellar multiple systems are underrepresented 
in our sample. However, since the sample size of our As- 
tralux survey is not yet statistically significant (^ 231 for 
95 % confidence level and an error margin of 5 %) we can not 
draw any conclusion for the whole exoplanet population. For 
a detailed statistical analyzis of the properties of all known 
stella r multiple systems harboring extrasolar planets, we re- 
fer to lRoU et ahl l|2011. submitted!) . 

We will continue our current monitoring campaign in order 
to determine the multiplicity status of all exoplanet host 
stars, which then will yield the true multiplicity rate of 
planet-bearing stars. This will also eventually allow us to 
draw conclusions about the frequency of planets in multi- 
ple stellar systems, as well as their properties compared to 
planets which reside around single stars. 
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Table 7. List of all stars with no additional stellar components detected. We give the numbers of frames observed, each frame with an 
exposure time of 29 .54 ms. We also give an upper limit for the minimum mass detectable at 0.5, 1, 2 and Sarcsec using the models by 
iBaraffe et ahl jiggsj ) and assuming an average age for each epoch as listed in Table [5] 



star 




RA 




DEC 


Epoch 


^ of frames 




[Mq] at 


0.5 






at l" 




: [Mq] 


at 2" 






r" 

at o 


HD 1461 


00 


18 


41 


,8 


-08 03 


10.8 


14/07/2010 


55000 


0, 


157 


± 


0. 


oil 


0. 


,099 


± 


0, 


,003 


0. 


,083 


± 


0, 


,001 


0, 


,082 


± 


0.001 


BD-1763 


00 


28 


34, 


,3 


-16 13 


34.8 


07/09/2009 


50000 


0, 


110 


± 


0, 


006 


0, 


,091 


± 


0, 


,003 


0, 


,083 


± 


0, 


,001 


0, 


,079 


± 


0.001 


Hat-P-19 


00 


38 


04, 


.0 


4-34 42 


41.6 


27/07/2011 


50000 


0, 


140 


± 


0, 


022 





109 


± 


0, 


.012 





,099 


± 


0, 


.008 


0, 


,098 


± 


0.007 


HD 4203 


00 


44 


41 


,2 


+20 26 


56.1 


11/07/2008 


80000 





176 


± 


0, 


020 





117 


± 


0, 


,008 





,091 


± 


0, 


.003 





,086 


± 


0.002 


HD5319 


00 


55 


01, 


.4 


+00 47 


22.4 


11/07/2008 


50000 


0, 


424 


± 


0, 


047 





221 


± 


0, 


.027 


0, 


140 


± 


0, 


.013 


0, 


128 


± 


0.011 


HIP 5158 


01 


06 


02, 


,0 


-22 27 


11.3 


14/01/2011 


50000 





112 


± 


0, 


008 





,089 


± 


0, 


.002 





,080 


± 


0, 


.002 





,079 


± 


0.001 


HD6718 


01 


07 


48, 


.6 


-08 14 


01.3 


14/01/2011 


50000 





211 


± 


0, 


021 





107 


± 


0, 


.005 


0, 


,085 


± 


0, 


.001 


0, 


,083 


± 


0.001 


HD 7924 


01 


21 


59, 


.1 


+76 42 


37.0 


27/07/2011 


50000 





100 


± 


0, 


003 





,082 


± 


0, 


.001 


0, 


,079 


± 


0, 


.001 


0, 


,078 


± 


0.001 


HD8673 


01 


26 


08, 


,7 


+34 34 


46.9 


14/01/2011 


50000 


0, 


,270 


± 


0, 


025 





113 


± 


0, 


.006 


0, 


,089 


± 


0, 


.002 





,088 


± 


0.001 


HD 9446 


01 


33 


20, 


,1 


+29 15 


54.5 


14/07/2010 


55000 


0. 


148 


± 


0, 


Oil 


0, 


,099 


± 


0, 


.003 


0, 


,086 


± 


0, 


.001 


0, 


,083 


± 


0.002 


HD 16232 


02 


37 


00, 


,5 


+24 38 


49.9 


23/02/2010 


20000 


0, 


,402 


± 


0, 


033 


0, 


,318 


± 


0, 


.034 


0, 


,208 


± 


0, 


.022 


0, 


121 


± 


0.007 


HD 16175 


02 


37 


01, 


,9 


+42 03 


45.4 


27/07/2011 


50000 





,236 


± 


0, 


024 





108 


± 


0, 


,005 





,093 


± 


0, 


.003 





,091 


± 


0.003 


HD 16760 


02 


42 


21 


.3 


+38 37 


07.2 


23/02/2010 


20000 


0, 


,289 


± 


0, 


047 





223 


± 


0, 


.036 





159 


± 


0, 


.019 


0, 


108 


± 


0.009 


HD 17092 


02 


46 


22, 


.1 


+49 39 


11.1 


14/07/2010 


55000 


0, 


441 


± 


0, 


181 


0, 


213 


± 


0, 


.122 


0, 


121 


± 


0, 


.043 


0, 


112 


± 


0.035 


HIP 12961 


02 


46 


42, 


.8 


-23 05 


11.8 


14/01/2011 


50000 


0, 


105 


± 


0, 


004 


0, 


,077 


± 


0, 


.001 


0, 


,074 


± 


0, 


.001 





,074 


± 


0.001 


HD 17156 


02 


49 


44, 


.4 


+ 71 45 


11.6 


23/02/2010 


20000 


0, 


,465 


± 


0, 


,038 


0, 


,376 


± 


0, 


.038 


0, 


243 


± 


0, 


.027 





136 


± 


0.010 


HD 22781 


03 


40 


49, 


.5 


+31 49 


34.6 


27/07/2011 


50000 


0, 


112 


± 


0, 


,005 


0, 


,086 


± 


0, 


.001 


0, 


,079 


± 


0, 


.001 





,078 


± 


0.001 


HD 28305^ 


04 


28 


36, 


.9 


+ 19 10 


49.5 


07/09/2009 


10000 


0. 


,862 


± 


0, 


,031 


0, 


,653 


± 


0, 


.028 


0, 


342 


± 


0, 


.036 


0, 


,325 


± 


0.034 


HD 32518 


05 


09 


36, 


.7 


+69 38 


21.8 


14/01/2011 


50000 


0. 


,686 


± 


0. 


,036 


0, 


,259 


± 


0, 


.030 


0, 


183 


± 


0, 


.020 


0, 


185 


± 


0.020 


HD 34445 


05 


17 


40, 


,9 


+07 21 


12.0 


14/01/2011 


50000 


0, 


219 


± 


0, 


,020 


0, 


,094 


± 


0, 


.003 


0, 


,089 


± 


0, 


.002 


0, 


,089 


± 


0.002 


HD 33564 


05 


22 


33, 


.5 


+ 79 13 


52.1 


07/09/2009 


50000 


0, 


,253 


± 


0, 


022 





148 


± 


0, 


,009 





097 


± 


0, 


.003 





,095 


± 


0.002 


HD 38801 


05 


47 


59, 


,1 


-08 19 


39.7 


14/01/2011 


50000 





,363 


± 


0, 


,065 





115 


± 


0, 


.012 





,097 


± 


0, 


.006 





,095 


± 


0.005 


HD 45652 


06 


29 


13, 


,2 


+ 10 56 


02.0 


16/01/2009 


22610 





165 


± 


0, 


,013 





118 


± 


0, 


.007 


0, 


,094 


± 


0, 


.003 





,087 


± 


0.001 


HD 45410 


06 


30 


47, 


,1 


+58 09 


45.4 


14/01/2011 


50000 


0, 


474 


± 


0, 


,034 


0, 


198 


± 


0, 


.019 


0, 


127 


± 


0, 


.008 


0, 


124 


± 


0.007 


HD 68988 


08 


18 


22 


,1 


+61 27 


38.5 


23/04/2008 


50000 


0, 


177 


± 


0, 


,015 


0, 


120 


± 


0, 


.007 


0, 


,092 


± 


0, 


.003 


0, 


,084 


± 


0.002 


HD 73534 


08 


39 


15 


.8 


+ 12 57 


37.3 


14/01/2011 


50000 


0, 


,287 


± 


0, 


041 


0, 


107 


± 


0, 


.007 


0, 


,093 


± 


0, 


.004 


0, 


,093 


± 


0.003 


HD 74156 


08 


42 


25, 


.1 


+04 34 


41.1 


16/01/2009 


16000 


0, 


,386 


± 


0, 


,040 


0, 


216 


± 


0, 


.023 


0, 


141 


± 


0, 


.011 





131 


± 


0.009 


HD 75898 


08 


53 


50, 


.8 


+33 03 


24.5 


23/04/2008 


50000 


0. 


222 


± 


0, 


,027 





130 


± 


0, 


.010 





,096 


± 


0, 


.003 





,091 


± 


0.003 


HD 81688 


09 


28 


39, 


.9 


+45 36 


05.3 


14/01/2011 


48149 


0, 


,647 


± 


0, 


,032 


0, 


291 


± 


0, 


.034 


0, 


216 


± 


0, 


.024 


0, 


216 


± 


0.024 


HD 87883 


10 


08 


43, 


,1 


+34 14 


32.1 


14/01/2011 


50000 





107 


± 


0, 


,004 





,084 


± 


0, 


.001 





,078 


± 


0, 


,001 





,078 


± 


0.001 


HD 90043 


10 


23 


28, 


.3 


-00 54 


08.0 


14/01/2011 


50000 


0, 


553 


± 


0, 


,032 


0, 


183 


± 


0, 


.017 


0, 


125 


± 


0, 


.008 


0, 


126 


± 


0.008 


HIP 57050 


11 


41 


44, 


,6 


+42 45 


07.1 


14/01/2011 


50000 




< 


0.072 




< 


0.072 




< 


0.072 




< 


0.072 


GJ436 


11 


42 


11 


,1 


+26 42 


23.7 


11/07/2008 


80000 


0, 


,080 


± 


0, 


,001 





,076 


± 


0, 


.001 


0, 


,073 


± 


0, 


.001 




< 


0.072 


HD 104985 


12 


05 


15 


.1 


+ 76 54 


20.6 


16/01/2009 


26480 


0, 


,682 


± 


0, 


,030 


0, 


,518 


± 


0, 


.034 


0, 


,253 


± 


0, 


.025 


0, 


,238 
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the companion candidate llMason et al.ll2009l'l at 0.237 arcsec separation could not be resolved 
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Figure 14. Average dynamic range in all observation epochs af- 
ter PSF subtraction. All detected stellar companions have been 
added. The typical standard deviation is 0.2 mag. 
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